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The usual idealized model for a cylindrical antenna consists of an extremely thin-walled 
tube of infinite conductivity with an infinitesimally narrow circumferential ^ap. Physically, 
one expects the antenna current at the feedpoint of such a model to be infinite. However, 
the singularity in feedpoint current is not detected by either iterative or Fourier series 
solutions of the integral equation for antenna ciuTent. These solutions are continuous 
on \z\ <h, where h is the half-length of the antenna. They are also in good agreement with 
each other and with experimental data. From a formal ])oint of view this amounts to saying 
that the conventional solutions of the integral eciuation are solutions on most, but not all, 
of the range j^l <h. Inside a small region, \z\<Zs, the correct solution to the integral equa- 
tion is peculiar to the idealized generator. Low-order iterative and P^ourier series solutions 
ignore these peculiarities and produce a smooth current distribution which can l)e used 
in predicting the behavior of practical structures. This paper is concerned with a de- 
tailed study of the theoretical current near the feed|M)int of the idealized model and \\w 
development of a definition of z^. 



1. Introduction 



The current distribution, It{z), for a finite cylin- 
drical antenna of length 2h and radius a is a solution 
of 
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R = [^a? ^xn^ <t>l2 + {z-z,y]K (3) 

Some writers approximate the kernel distance with 

R^=W+{z-z,y]K (4) 

With the correct kernel distance the integral equation 
corresponds exactly to a model consisting of a thin- 
walled tube excited by a hypothetical generator of 
voltage Fsuch that the electric field in a narrow gap 
is given by 

E,,,=-Vb{z). (5) 

Idealization of the generator region leads to a 
fundamental difficulty. A finite voltage across an 
infinitesimal gap implies that I{z) is infinite at z=0. 
However, it is possible to manage the solution of (1) 
in such a way as to obtain a useful finite result. In 
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the King-Middleton [King, 1956] and Hallen [1938] 
theories the integral equation is solved by an iter- 
ative process in which the starting approximation is 
a continuous function. After only one or two 
iterations the output function is still continuous. 
It can then be argued that these solutions are valid 
on \z\<h except for a small region \z\<Zs where they 
differ markedly from the mathematically correct 
solution. The continuous ''solution'^ evaluated at 
z=0 gives a finite admittance, Ya, which is presumed 
to be characteristic of the antenna viewed as a 
circuit element. There is still the problem of 
establishing correspondence with experiment. 
Measiu^ed antenna admittances depend on Ya and 
tlie actual terminal zone structiu'e and transmission 
line used in the experiment. Experimental pro- 
cedures and theoretical corrections can be designed 
to extricate Ya from the experimental data [King, 
1956]. Combined theoretical and experimental ex- 
perience have definitely established the utility of 
approximating the correct solution of (1) by a 
continuous function. 

Iterative solutions of (1) become extremely tedious 
beyond second or third order*. Some of the integrals 
arising in the theory are performed with the approxi- 
mate kernel distance given by (4). This approxima- 
tion is suspect for thick structures. These objections 
are overcome in a theory developed by Duncan and 
Hinchey [I960]. This theory provides a method for 
obtaining the coefficients of a Fourier expansion of 
the current distribution to any reasonable degree of 
accuracy. The solution has been carried to 25th 
order. One interesting result is that no sign of poor 
behavior of I{z) near z=0 develops even at such high 
order. One does not expect a Fourier series to show 
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fine details in a region of width Zs until the order of 
solution is high enough for space frequencies of order 
2^^ to appear. The fundamental period in finite 
antenna theory is about a wavelength. Thus, our 
calculations provide evidence that 

0,<X/25for A/a>60, (6) 

the smallest h/a ratio studied. 

According to Wu and King [1959], Zs may be much 
smaller than the estimate provided by (6). It can 
be shown that 



K^y- 



. ikaV 1 , 
-J —y — in kz 



(7) 



for small z. If one is allowed to guess that this func- 
tion is simply not applicable after it decreases to a 
value comparable to the amplitudes of well-known 
continuous solutions one obtains 
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Although (7) is certainly correct for small enough z, 
the method of obtaining (8) is somewhat arbitrary. 
A definition of Zg based on the point at which a simple 
logarithmic function exceeds the nominal value of 
observed current fails to take into account the man- 
ner in which the current makes the transition from 
moderate to logarithmic behavior. The definition 
of Zs to be introduced here is a rough measure of the 
location of the transition region in which the graph 
of the current changes from moderate to large slope. 
On this basis (8) turns out to be much too small and 
too sensitive to antenna diameter to be useful. 

Obviously, the distinction between thin and thick 
cylinders will always be somewhat arbitrary. Such 
distinctions depend on methods of analysis and de- 
sired numerical tolerances. In this paper a thin 
antenna is one in which a moderately smooth func- 
tion is a good approximation to I{z) except within a 
suitably small fraction of a wavelength near 0=0. 
The remainder of the paper is devoted to the quan- 
titative development of this general idea. 

In order to facilitate the analysis we allow the 
antenna length to become infinite. The singularity 
in I{z) is associated entirely with the idealized gen- 
erator. Waves reflected from the ends of a finite 
antenna can do no more than produce a standing 
wave of twice the amplitude of a smooth approxima- 
tion to I{z) for the infinite antenna. 

A byproduct of this study is the development of a 
continuous approximation to I{z) for the infinite an- 
tenna when ka is small enough to justify such a cur- 
rent distribution. Admittances defined from the 
approximate solutions for various values of ka are to 
be interpreted similarly to those of finite antenna 
theory. Admittance values for infinite antennas are 
of more than academic interest. Altshuler [1961] 
recently reported on a traveling-wave linear antenna. 
Lumped resistors are inserted in this type of antenna 
at \z\ = h — \l^. With appropriate values of resist- 



ance the current distribution in the central region, 
0<^|s|<^A — X/4, consists primaril}^ of running waves 
traveling outward from z^O. It is thus possible to 
appl}^ the theor}' of the infinite antenna to a realizable 
structure. 

2. Formal Theory 

The right-hand side of (1) represents the surface 
value of vector potential. When the antenna length 
becomes infinite the term C cos kz which represents 
standing waves is no longer needed. With V=l 
the integral equation becomes 






lt{2,)g{z- 






Hallen [1956] has reported the solution of (9). 
the notation of this paper. 
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It{z) may be broken into two parts 
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The first of these expressions gives the current 
distribution on the outer surface of the tubular 
model. Current on the inner surface is described 
by (14). Only the exterior solution is of interest in 
antenna theor3^ Incidentally, (13) is in agreement 
with the work of Silver and Saunders [1950] who 
solved the exterior boundary value problem of an 
infinite C3dinder with a window of arbitrary size and 
shape in which the electric field is known. Equation 
(13) follows from their work if the window is allowed 
to become a narrow circumferential gap and the 
usual boundary condition relating surface current 
density and magnetic fiekl is imposed. 

The task now at hand is to integrate (13) in such 
a way that I{z) may be examined arbitrarih" close 
to ^=0. It is helpful to break the integrand into 
real and imaginary parts and write 

/(3)= r [R{a)+jW(oi)] COS azda, (15) 

I(z)--Mz)+ji{z), (16) 



or 
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where 0^=¥—a^. We have chosen the contour of 
integration along the real axis with an upward in- 
dentation at a=k. The phase of fi is defined b}^ 

phfi^^O, a<ik; php= -t/2, ayk. (11) 

With the help of the identity 

J,{x)Y,{:x)-J,{x)Y,{x)=2lxx, 



where 

II {a) 



ik 



A'(«)=0, 
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a<Ck: 
a>k; 
a<k; 

tyk. 



(17) 



Thus, the real part of I{z) can be determined by an 
integration only on the finite range 0<a<^. The 
singularity in B{a) at a=k is integrable so that 
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The combination of Bessel functions seen in (18) 
occurs so frequently in what follows that it is con- 
venient to define 



B{x)==Jl{x)+Yl{x). 



(19) 



Further reduction is accomplished by the change of 
variable, a=k cos 6. We also let kz^l; and ka=K. 
Then, 

cos (f COS B) 



'■^> xZoJo si, 



sill eD{K sin 6») 
Now intioduco the well-known identity 



de. 



(20) 



cos (f cos 0)= Z; ««(-!) V2„(r) cos (2ne); (21) 

en=l for n = Oy and e^=2 for ?i>0. 

In terms of coeflicients w^hich can be obtained b^^ 
numerical or approximate integration. 



with 
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Further discussion of the i?2n <'ind a table of values 
can be found in the section entitled Numerical 
Analysis. It is sufficient at this point to note that 
only a few terms of (22) are required to provide an 
accurate representation of r(f) on the range of 
0<f<27r. Since 1(2) is an even function, this 
corresponds to a full wavelength on either side of 
the feed point at 2=0. 

The determination of i(z) is more difficult. By 
using the small argument formulas for the Bessel 
functions and modified Bess(^l functions which occur 
in (17) it can be verified that \V(a) contains non- 



integrable singularities of opposite sign on either side 
oi a=k. A direct attack would involve evaluating 



i(z) =-r- W(a) cos azda, 
Jo 



(24) 



where a standard notation has been used to denote 
the principal value of the integral. Defining i(z) by 
means of the principal value of an ini})roper integral 
is not hard to justify. The difHculty could have 
been avoided in the first place by detouring well 
away from a = k in defining the contour. Even 
with this admitted, a serious (hfiiculty remains. 
W(a) is a complicated function and the integral must 
be performed numerically over a long range in a 
before Ki and Kq are accurately represented by their 
asymptotic behavior. To do this and maintain high 
accuracy for small z is (Ufficult, even with the aid of 
automatic com])uter facilities. 

Fortunately, the difficulties inherent in using (24) 
can be avoided entirely. R(a) can be used to deter- 
mine both r(z) and i{z). We rewrite (13) as 

/^ 00 

I{z)=\ F(a) C0-^aZ(Ja. (25) 

An obvious substitution gives 

I(z) = -j(^ F(a)^maZ(la+( F{a)eJ^'da. (26) 
Jo Jo 

Denoting the ijuaginary |)arl of the first integral by 

^i(^), 

%,{2) = - r/t(«) sin azda. (27) 

A series form similar to (22) can l)r written using 

sin (f cos 



2 2: (-l)"J2„+,(f) cos (2n^-l)e. (28) 

77 = 



One obtains 
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OS {2n+l)d 
eD(K sin d) 



dd' (29) 



A table of /i and Iz for various values of K can be 
found in the next section. 
There remains 



'i;2(^)=Im 



F(a)e^'''~da. 



(30) 



We now choose a closed contour consisting of a line 
from the origin along the positive real axis with an 
upward indentation at a=k, the infinite quarter 
circle from the real axis to the positive imaginary 
axis, and a return to the origin via the imaginary 
axis. If 2 is positive the integral along the large 
quarter circle will vanish. There are no poles within 
the contour. Then 
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^2(2)=Im 



F{a)e^'^hla. 



(31) 



The substitutions a=jr], 7=V^^+V reduce (31) to 

ttZo Jo Y^iyd) 

Even ^2(.3') is determined by R{a). With the under- 
standing that R(ja) is defined on an infinite range, 
whereas E{a) — for a^k, we can summarize in 
terms of cosine, sine, and Laplace transforms with z 
as the transform parameter. Equations (18), (27), 
and (32) can be collected into 



Re 7(3)=.^/? (a), 

Im/(^) = -.^/?(a)+^'7?(ia). 



(33) 



In this sense the real and imaginary parts of I{2) are 
codetermined. 

For computational pui'poses (29) and (32) are a 
considerable improvement over (24). However, (32) 
can be improved further. The first step is the triv- 
ial one of obtaining an integrand in terms of non- 
dimensional quantities. If we let rj = kv and ka=K, 



i2(r)= 
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(l+v')DiK^ll+v') 
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With the substitution u=K^^l+v'^ (34) is trans- 
formed to 

i2{^)=^ I , ;^;r^^x ^ :~du. (35) 
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The advantage of (35) over (34) is apparent from 
the asymptotic behavior of [uD{u)]~\ At ^2^=10 
this function is approximated b}^ 7r/2 with an accuracy 
of 0.12 percent. Thus we need be concerned with 
numerical integration only on the range K<u<10 
no matter what the value of K. By contrast, the 
integrand of (34) cannot be replaced by its asymp- 
totic behavior until vc^lO/K. For large K (35) can 
be integrated directly [Erdelyi et al., 1953] to give 



^.(f)~^[H„(f)-^„(f)], 



(36) 



where Ho(f) is the Struve function of zeroth order. 
A difterent strategy is to be employed for small K. 
To eliminate unnecessary repetition of constants 
during the derivation we define K times the integral 
in (35) as /(f). We start with the identity 
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The first of these is integrable. In the third integral 

we replace {uD{u)]~^ by 7r/2 and neglect Kin^lu^—K^. 
The second term is defined as a special function, 
qxi^), to be determined and tabulated by numerical 
integration. We now have 



juy- 



-2D{K) 



[Ho(f)-F„(f)]-fe(f) 



Finally, 



+\j^-^~\Ei{-miK). (38) 
i (f ) =^ \j (f ) +'l5 /2„+. J2.+l(f )1 (39) 



unless K>10, in which case (4/7rZo)/(f) is replaced 
by (36). 

3. Numerical Analysis 

The coefficients i?2n a.nd /2n+i have in common a 
factor of the type 



"■=1 



cos rW 



sin eD {K sin 6) 



de. 



(40) 



Near ^=0 the integrand of (40) becomes approxi- 
mately 

[e[\ + {4iir')\nHyKei2)]]-\ (41) 

where In 7=0.57722. It is convenient to break the 
integration into two parts. In the first part 0<^^<;€ 
and the integrand is replaced by the small argument 
formula (41). The singularity is integrable, a 
convenient substitution being {21 it) ln(yKd/2)=y. 
The result of integrating formalh^ in the small d 
region is 



^g+arctan01n^)] 



Numerical integration was used in the remaining 
interval, €<^<7r/2, and the result added to (42). 
The parameter e was always chosen so that the 
small argument formula for Fo(^sin 6) was accurate 
to four significant figures. A short table of i?2?^ and 
l2n+i coefficients, adequate for the computation of 
/(f) on the range 0<f<l, is presented in table 1. 







Table 1 




K 


7?o 


/i 


R2 


Iz 


0.01 


0. 5683 


-1.009 


-0. 7759 


0. 6707 


.02 


.6678 


-1. 168 


-. 8631 


.7283 


.04 


.8054 


-1.383 


-.9694 


.7912 


.08 


1.005 


-1.686 


-1.099 


.8550 


.15 


1.277 


-2. 086 


-1.244 


.9035 



Tables of qxi^), to be omitted here for the sake of 
brevit}^, are being prepared for separate publication. 
However, a few of the technical details and a quali- 
tative description of the results are of general 
interest. Recall that g/v:(f) is defined by 
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i^lK) -ylu-^-K^ 



-du. (43) 



K' 



The significant features of the design of gA'(f) are: 

(a) For i^=10, g^CD-O. If i^>10, qM) is not 
needed at all since the contribution to ^(f) from the 
imaginary axis integration is adequately represented 
by (36)/ 

(b) The integrand defining q^Ki^) is everywhere 
finite on a finite range of integration. This feature 
simplifies, to a great extent, the problems of accuracy 
control and computer programing. 

A plot of a typical qK(^), figure 1, shows that this 
function is appreciable near the origin, concave 
upward, and highly damped. 

The approximations introduced in the third term 
of (37) are not serious even for large K. Corrections 
can be derived by expanding the integrand in a 
serie s an d integrating term by term. Setting K=0 
in -yju^—K^ simply provides the first term of the 
series. Correction terms turn out to be of the form 
t''Ei{—\Q^IK). Near f-0 the correction is small. 
Near f=l the correction amounts to about 50 per- 
cent of the last term of (38) for K=5. The error 
decreases for smaller K. However, for f near ^=1 
the contribution to the total current from the integral 
in question is small. The net effect is a small error 
in a large term near the origin passing over to a 
Lirge error in a small term as f is increased. For- 
tuitously, the correction terms may be dropped 
entirely even when -y^^u^—K^^^u is a poor ap])roxi- 
mation. 

The singularity in i(^) is expressed explicitly by 
ro(f) and Ei{-m/K). The graph of i(r) becomes 
infinite at ^=0, decreases rapidly in a short interval, 
and then changes in a small transition interval to a 
moderate slope. As f increases ^^(f) decreases and 
the Bessel functions in i»*i(f) begin to dominate and 
establish the wave character of i(f). Plots of i(f) 
were made for various values of K. A typical plot 
for K=O.Ol is shown in figure 2. The singularity 
in the current distribution does not appear in the 
theory of antennas of finite length as developed to 
date. One would like an analogous current distri- 
bution in the theory of the infinite antenna. We 
have followed the intuitive procedure of extrap- 
olating from the slowly varying portion of ^(f) into 
the region of small f. The resulting graph with the 
rapidly varying portion of ^(f) removed is designated 
as ^s(f)• The range of the singularit}^ is arbitrarily 
defined bv 



i{^s) = is(0); z,l\=tJ2ir. 



(44) 



For small K, Zs is a small fraction of a wavelength. 
Values of 2s estimated from graphs of i(f) for various 
values of K are presented in table 2. The first line 
for 7^=0.001 was estimated from the behavior of the 
remaining data. As i^is increased the solution turns 
yery gi-achially into the logarithmic singularit}^ 
Figure 3 exhibits two extremes with plots of ^(f) for 




Figure 1. Graph of r/K(f) versus f for X=0.01 
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Figure 2. Expanded graph of i{^) versus ^ for i!C= 0.01 
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Table 2 



t — ■— — ■- 
K 


U 


Zsl\ 


0.001 


0.006 


0. 00095 


.01 


.008 


.0013 


.02 


.010 


.0016 


.04 


.015 


.0024 


.08 


.020 


.0032 


.15 


.024 


.0038 




Figure 3. Graphs of i{^) for K=0.01 and K=10 



ir=0.01 and i!L=10. The singularity in the graph 
for ir=0.01 is not evident unless the scale of the 
abscissa is expanded considerably. 

It is at once apparent that (8) is not an adequate 
estimate of the range of the singularity. The tran- 
sition from moderate to large slope in the graph of 
'i(f) cannot be accounted for entirely by a simple 
logarithmic function. The qualitative conclusion 
that iterative and Fourier series solutions of feasible 
order will not detect the singularity remains intact. 



4. Modification for Nonzero Gap Width 

According to Silver and Saunders [1950] a gap in 
the region —W<z<W is accounted for by a factor 
(sin at^\)/a:TT^in the integrand of (18). If this factor 
is included in the present work it will have very 
little effect on the functions r(f) and ii(f). If the 
gap width is a small fraction of a wavelength (sin aW)l 
aW will be nearly unity throughout the range of the 
integrals which define these functions. The appro- 
priate modification turns out to be the inclusion of a 
factor (sinh vw)/vw in the integrand of (34), giving 






y-vt 



sinh vw 



(l+v')D{K^Jl+v^) ^'^ 



dv, (45) 



where w=kW, Denote the unmodified integrand 
by F{v) for convenience. Then 

ic-)(f)=A r"^[6-^f-"^^-6-^f+"^^]rf?;, (46) 



wv 



^2 



TTZ/q Jo 



'^dv, 



(47) 



where appropriate superscripts have been used to 
identify the separate cases. It is easy to show that 



1 ff+«' 



>(f)^f. 



(48) 



Thus, the required modification is a simple average 
over the evanescent part of the current distribution 
for a 5-gap. Results for the two cases will be appre- 
ciably different only for small f . It is more difficult 
to specify how this solution is to be matched to that 
for an appropriate transmission line and terminal 
zone structure. In thin antenna theor}^ this diffi- 
culty is resolved by treating the antenna and termi- 
nal zone structure separately by suitable approxi- 
mations and combining the results in an equivalent 
circuit. Clean separation of these two parts of the 
problem is not possible if the rapidly var^ ing portion 
of ^2^^ (r) extends over too large a region. 

The simplicity of (48) can be attributed to the 
assumption of constant electric field in the gap. 
Wait [1959] has pointed out that a more realistic as- 
sumption leads to a factor of Jo(aw) in the integrand 
of (18) instead of (sin aw)/aw. 

5. Impedance 

The conductance of the cylindrical antenna is eas- 
ily defined for all values of K. It is given simply by 



G^riO)-- 



_4 



'oJo si 



dd 



sin e D(K sin 6) 



(49) 



The integral is the definition of the coefficient i?o 
defined earlier. If K is large the function D(K sin^) 
takes on its asymptotic value of (tK sin d)/2 through- 
out most of the range of integration. Then the 
conductance-circumference ratio is 
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g=G/2Tac^7r/Z^X 



(50) 



This is, as it should be, one-half the coiKhictaiicc^ per 
unit len<2:th of a slit in an infinite plane [Jordan, 
1950]. Only one side of the cylindrical surface is 
counted because of the original splitting of the total 
current into antenna current and tube current. At 
K='{0 (50) gives a value of 0.0835 mhos compared 
to a correct value of 0.0876 mhos. Cylinders with 
K> 1 are essentially infinite in diameter for the 
purpose of admittance calculations. 

Susceptance is defined from the smoothed version 
of i(r) by 



B= 



vz; 



^.(0). 



(51) 



This represents an attempt to estimate, on a plau- 
sible basis, the value of susceptance observed in 
experimental work. Susceptances found in this man- 
ner are analogous to those calculated by the King- 
Middleton theory or Fourier series solutions for the 
anteiuia of finite length. As K is increased the 
introduction of a smoothed version of i{^) becomes 
less tenable. We have arl)itrarily selected K=0.15 
as that value above which the concepts of thin an- 
tenna theory are not to be api)lied. This is simply 
a matter of judgment based on the Ix^havior of curves 
of '?(r) versus f for various values of K. 

For K>10 a gap width, 2W, can be us(mI in con- 
junction with (36) and (48) to define 



B^ 



;{w) 



(w)- 



2wJo 



\t)<lt. 



(52) 



Only the Neumann function in (36) will contribute 
appreciably to (52). Passing to the small argnnuMit 
formula for Fo(f) JViul carrying out tlu^ indicated 
operations gives 

J = 5/2.a=.^(ln^+l) (53) 

This is the susceptance per unit length for a slit in 
an infinite plane when the excitation is defined as a 
constant electric field in the slit. The slight differ- 
ence between (53) and the corresponding formula in 
Jordan [1950] can be attributed to the differences in 
assuuKHl excitation. 

Table 3 gives the calculated admittances and im- 
pedances for small K. 

Table 3 





Admittance 


Impedance 


K 


GXW 


J5X103 


R 


X 


0.01 


1.92 


0.761 


450 


-178 


.02 


2.26 


1.08 


360 


-172 


.04 


2.72 


1.55 


278 


-158 


.08 


3.40 


2.43 


195 


-140 


.15 


4.32 


4.39 


114 


-116 



Altshuler [1961] recently reported the measured 
impedance of a traveling wave linear antenna with 
7^=0.04 as (320-j 110) ohms. The transmission 
line in his experiment was a coaxial line with a I'atio 
of outer to inner diameters of the conductors of 8.13. 
The predicted value based on the present theory and 
appropriate terminal zone coi'rections is (346-^* 86) 
ohms. 

6. Conclusions 

The following classification of cylindrical antenna 
problems according to values of K, the antenna 
circumference in wavelengths, is suggested by the 
present theory and the literature. For K <0.1 a 
smoothed version of the exact solution of llallen's 
integral equation can be used together with appi'o- 
priate equivalent networks for realizable terminal 
zon(^s to predict the behavior of physical antennas. 
Smoothing of the solution involves only a small 
region near the lecMlpoinl of the idealized mo(l(4 to 
which the integral ecjnation apj)lies. The size of this 
region is of the order of a thousand 111 of a wave- 
length. This is much larger than had been previ- 
ously supposed but still suffici(Mitly s.mall to justif}^ 
the basic approach. Appi'oximations to the kernel 
of the integral equation can be safely used for 
/^<0.02. In particular, the King-Middleton theory 
is in agreement with high-order Fouiier-series 
solutions using the exact kernel. For 0.1<^7C<;i0 
the antenna and terminal zone problems are not 
easily sepai'ated. A th(H)ry to j)rovide botli im])ed- 
ances and radiation charactei-istics for tliese (ham- 
eters has not yet been developtMl. An adequate 
theory will have to involve a unified treatment of the 
terminal zone and a gap of nonzero width. Above 
iiL>10 the impedance properties of a long cylindrical 
antenna are well ap])roximate(l by those for a slit in 
an infinite plane. 

The theory and lesults presented here ])rovide 
information about the singularity in theoi'etical 
antenna current. The calculated impechmces should 
prove useful in the design of traveling- wave linear 
antennas. 



The author thanks K. S. Kunz and C. W. Harri- 
son, Jr. for continued encouragement and helpful 
discussions. J. F. Stanley and S. E. Dickson 
performed most of the calculations. 
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